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Finite volume effects in chiral perturbation theory 

Gilbcrto Colangclo a 
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Sidlerstrasse 5, 3012 Bern, Switzerland 

There has recently been an intense activity in the study of finite volume effects by means of chiral perturbation 
theory. In this contribution I review recent work in this field both for the e- (M^L < 1) and the p-regime 
(M n L ^> 1). For the latter I emphasize the importance of going beyond leading order calculations in chiral 
perturbation theory and the usefulness of asymptotic formulae a la Liischer used in combination with CHPT. 
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1. Introduction 

Every lattice calculation is done in finite vol- 
ume and at finite lattice spacing, as well as at 
finite (and at present typically substantial) quark 
masses. Before being able to make a meaningful 
comparison to experimentally measured quanti- 
ties, one has to make three extrapolations - and 
since all of them are numerically quite expensive, 
any analytical method that may help in this re- 
spect is more than welcome. Chiral perturbation 
theory (CHPT) provides the proper framework 
for calculating analytically the dependence on all 
three extrapolation parameters. An overview of 
the present status of the analytical calculations 
related to the extrapolation to the continuum and 
to the chiral limit, and in particular of the inter- 
play between the two limits, has been given by 
Oliver Bar £Q. Here I will discuss the extrapo- 
lation to the infinite volume limit, and give an 
overview of the present status of analytical calcu- 
lations in the framework of CHPT. 

If a physical system is enclosed in a finite box, 
and periodic boundary conditions are imposed, 
the space of the possible three-momenta becomes 
discrete: 

p=^n, n = (ni,n 2 ,n 3 ) (1) 

with ni,2.3 S Z. If the box is large enough, this 
discrete space will almost look like a continuum 
one: and since one is interested in simulating 
the soft, nonperturbative dynamics of QCD, one 
would like to have the region of soft momenta to 



look like a continuum one. If we take 4ttF v as the 
chiral symmetry breaking scale, which separates 
soft from hard momenta, we obtain the following 
quantitative condition on the volume 



2tt 

— < 4irF„ 



X »^~ lfm ( 2 ) 



in order to have a large number and not only a 
few discrete momenta in the soft region. 

If the condition J2J is satisfied one can then 
use CHPT to study analytically the behaviour 
of the system at low momenta, and in particular 
the explicit dependence of physical observables 
on the volume. How to extend the CHPT frame- 
work to the finite volume case has been discussed 
by Gasser and Leutwyler P|: in finite volume 
CHPT becomes a systematic expansion in both 
the quark masses and the inverse box size. In 
infinite volume there are relations among the co- 
efficients of this expansion for different quantities 
which follow from the chiral symmetry of QCD: 
these relations go traditionally under the name 
of "low energy theorems" and CHPT is a conve- 
nient tool to derive them systematically. In finite 
volume chiral symmetry implies relations among 
the coefficients of the expansion in 1/L of differ- 
ent observables as well as relations among these 
coefficients and those of the expansion in quark 
masses. CHPT allows one to obtain these "large 
volume theorems" in a systematic way. 

The product F^L is not the only relevant pa- 
rameter. An important role is also played by the 
relative size of M v and 1/L. If M^L < 1 finite 
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volume effects will be similar to those in the chiral 
limit - as is well known, no spontaneous symme- 
try breaking can take place in a finite volume, and 
one will see a deformation of the vacuum state due 
to finite volume effects. If M n L 3> 1, on the other 
hand, the effect of the explicit symmetry breaking 
on the low energy behaviour of the system will be 
more important than that due to the finite vol- 
ume. In the framework of CHPT, the difference 
between the two situations is given by the impor- 
tance that the Goldstone-boson zero modes have 
in the evaluation of the path integral. In the lat- 
ter case, like in infinite volume, the contribution 
of the zero modes can be neglected, whereas in 
the first case is important and must be explicitly 
evaluated |2I3| . The first situation is denoted as 
"e-regime" and the second one "p-regime", and 
the respective counting schemes are as follows: 

e-regime: M r ~ — r ~ (3(e 2 ) 

p-regime: M n ~ — ~ 0(p) . (3) 

In the first case the Compton wavelength of the 
pion is larger than the box size, and the pion does 
not have enough space to propagate (all other 
non-Goldstone-like particles however do, because 
their mass is of the order of A ~ 4:%F W , and we 
assume that the condition (0 is satisfied). In 
the latter case, a pion fits comfortably well in- 
side a box and, being the lightest particle, is the 
only one that has a nonncgligible probability to 
propagate until the box boundaries, and there- 
fore feel their presence. The different behaviour 
of the physical system in the two cases is well il- 
lustrated in Fig. n where the time dependence of 
the correlator of two axial charges 

(Qt(t)Qt(0))^S lk L 3 T A (t) (4) 

is plotted for different values of M W L. In infinite 
volume the correlator drops exponentially with a 
rate which is proportional to M^: 

T 1 T°"(t)=M„F*/2e- M - t [l + 0(p 4 )] . (5) 

In the p-regime one expects a behaviour of the 
correlator which is qualitatively similar to the one 
in infinite volume - as Fig.[I]rndeed shows. In the 
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Figure 1. Time-dependence of the axial-charge 
correlator on for different values of M^L. Note 
that for M^L = 2 the calculation in the e- and p- 
regime give perfectly overlapping results. Figure 
taken from |3] 



e-regime, on the other hand, the qualitative be- 
haviour is completely different from the infinite- 
volume case - the curves in Fig. ^ corresponding 
to small values of M V L can be well described by 
a polynomial of low degree in t [3] . 

2. Finite volume effects in the e regime 

One of the goals of lattice calculations is to de- 
termine the low-energy constants of the CHPT 
Lagrangian. These constants are by definition in- 
dependent of the light quark masses and a reli- 
able determination should therefore be made with 
quark masses as small as possible, which is noto- 
riously expensive. Even when very small masses 
will become accessible 1 working close to the chiral 
limit and in the p-regime will require enormous 
volumes and will likely remain for a long time pro- 
hibitively expensive. Approaching the chiral limit 
on the lattice implies working in the e-regime. As 
discussed above, observable quantities may look 
very different in the e-regime from the infinite- 
volume case. The strategy of the calculation in 

1 In the quenched approximation very small quark masses 
have already been simulated, see e.g. 0]. 
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this case must start from an analytical calcula- 
tion in CHPT to identify the proper observable 
that will allow an extraction of the low-energy 
constant of interest. The first analytical calcu- 
lations in this direction have been made in |3I5| . 
The first numerical calculations which exploited 
this strategy were aimed at the extraction of the 
quark condensate in the chiral limit |fil4| . 

More recently, two different groups engaged in 
the calculation of charge-charge correlators of the 
type defined in Eq. (J2J in order to extract the 
pion decay constant in the chiral limit. Results 
have been published during last year |7I8| and 
show a very clean determination of this chiral- 
limit observable. In fact one of the two groups has 
performed the calculation of the decay constant 
also in the p-regime and has shown that after the 
extrapolation to the chiral limit the result is in 
agreement with the direct determination in the 
e-regime: 

F ~ 130 McV 
F = (102±4MeV) e , (104 ± 2 MeV) p 

Although the calculations have been made in 
the quenched approximation, the positive results 
prove the feasibility of the method and have 
opened the way for a more ambitious aim of the 
program: the calculation of the low-energy con- 
stant of the weak chiral Lagrangian for nonlep- 
tonic if-decays. The first results in this direction 
have been presented at this conference and 
have also recently appeared [TU) . 

Unfortunately, these groups have also found 
out that working in the e-regime is technically 
challenging, and to complete their program had 
to overcome unexpected, highly nontrivial diffi- 
culties of numerical nature |7I8I9| . It will be in- 
teresting to see the application of this strategy 
with dynamical fermions - which is certainly a 
highly nontrivial task. 

Two papers have recently investigated this 
regime for the one-nucleon sector and promise in- 
teresting developments |11| . 

3. Finite volume effects in the p regime 

Gasser and Leutwyler f2| have shown that, if 
one works in the p-regime, in a large isotropic 



box with periodic boundary conditions, the only 
change in the calculation is that the propagator 
of the pseudo Goldstone bosons has to be made 
periodic: 

G L {x,t) = ^G OQ {x + nL,t) . (6) 

n 

The local vertices, as given by the CHPT 
Lagrangian, remain unchanged and the power 
counting for loop diagrams works like in infinite 
volume: every loop generates an additional 0(p 2 ) 
correction. The only difference is that in finite 
volume p ~ M ~ 1/L, and a one-loop calcula- 
tion gives both the leading quark-mass and finite 
volume corrections. For example, the one-loop 



calculation of M n 


and Fn gives 0: 


AUiL) 


= M v 


l + ^i(A) + 0(£ 2 ) 


FAL) 


= F„ 


Nt 

1- ^Si(A) + 0(£ 2 ) 


with A 


= M W L, £ 


= (Mrr 1 AttF-k ) 2 and 


.9i (A) 


= Elj 


f dze-i-^ 2 






k=l v feA 



is the tadpole correction in finite volume (m(fe) 
is the multiplicity of a vector fi with n 2 = k). 
The result shows that finite volume corrections 
for masses and decay constants are suppressed ex- 
ponentially, by increasing powers of exp(— M^L). 
The infinite sum appearing in Eq. |JSJ reflects the 
infinite sum in the definition of a fully periodic 
propagator, Eq. 10, but since M^L ^> 1 only the 
first few will be numerically relevant. 

A different approach which relies on an expan- 
sion in exponentials and not on CHPT is due to 
Luscher |12| - this expansion, however, has only 
been devised for masses and, more recently, for 
decay constants ^31 > an d is therefore less general 
than the chiral expansion. A detailed comparison 
of the two approaches will be given later. 

3.1. Overview of recent results 

The examples just discussed show that, in the 
p-regimc finite-volume effects are typically small 
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and relevant for very precise lattice calculations. 
Indeed only recently there has been an intense 
activity in analytical calculations of these effects 
despite the fact that the theoretical framework 
had been established around twenty years ago [2] 
12 , and extended to the quenched approximation 
not much later I14I15I. 



During last year there have been analyses of 
finite volume effects for: two-pion state energies 
[T5]. the pion mass J7|, Fk and Bk [TS]. the 
nucleon mass |19I20| . the nucleon magnetic mo- 
ments and axial coupling constant as well as the 
mass |2I], f B and B B [22] and F„ Many 
of these results have been also presented at this 
conference [23124125] together with new ones con- 
cerning M,r , F n and the charge radius [21] ■ Most 
of these calculations have been made in the frame- 
work of CHPT and concerned one-loop finite- 
volume effects and not using Luscher formula. 
A comparison to the Luscher formula has been 
made in the case of the proton mass ^U]i an d 
a discrepancy has been found with the explicit 
formula given by Luscher in his Cargese lecture 
notes [22] • It turns out that the formula in these 
lecture notes was not correctly written and that 
the correction has numerically important effects 

(cf. HHEOI). 

In what follows I will discuss in some more de- 
tail why it is important to go beyond the leading- 
order calculation of these effects if one wants to 
have a good control over them, and how this can 
be achieved by a combined use of CHPT and the 
asymptotic formulae a la Luscher. This is the 
approach advocated in |17ll8j . 

3.2. Hadron masses in finite volume: 
CHPT vs. Luscher Formula 

Above we have briefly discussed what rules one 
has to follow if one wants to calculate finite vol- 
ume corrections to a given quantity in CHPT, and 
shown the result for the pion mass. The analo- 
gous result for the kaon, the proton or any other 
hadron mass will of course be different. Luscher, 
however, has shown that the leading exponen- 
tial correction to the mass of a given hadron 
can be expressed in a completely universal way 
through an integral over the scattering amplitude 
of the hadron in question and the pion. The lat- 



ter enters here because it is the lightest hadron: 
the leading exponential correction to any hadron 
mass is of the form exp(— M W L), and how large 
the correction is depends on how strong is the in- 
teraction of the hadron in question to the pion. 
For example the Luscher formula for the pion 
mass reads 



dy F(iy) e -V M Z+v 2L + 0(e~ ML ) 



where F{v) is the physical (for v real - what en- 
ters the integral is the analytic continuation) for- 
ward 7T7T scattering amplitude and M > \plM. v . 
Note that the formula does not rely in any way 
on the chiral expansion, and indeed is valid also 
in cases where there is no spontaneous symmetry 
breaking and the corresponding Goldstone bosons 
- in such cases the role of the pion has to be over- 
taken by the particle which happens to be the 
lightest one. 

The formula can be used for a numerical anal- 
ysis only if one has a representation for the scat- 
tering amplitude in the integrand which lends it- 
self to a numerical evaluation. For hadrons one 
could in principle take directly the measured scat- 
tering amplitude (and do the necessary analytic 
continuation) to evaluate the integral, but then 
one would be restricted to physical values of the 
quark masses: in order to use the formula for 
arbitrary values of the quark masses one must 
therefore rely on the chiral representation for the 
scattering amplitude. Inserting the leading order 
chiral representation for the tttt scattering ampli- 
tude: 



M 2 



(10) 



in Eq. © one should obtain the same expression 
for the leading exponential term in Eq. 0. As 
was verified long ago |2], this is indeed the case. 

The chiral expansion for the scattering ampli- 
tude is now known up to and including order p 6 
[2"5]. and it is easy to insert the corresponding 
analytic expressions in Eq. © and evaluate the 
integral. This has been done in J7|, an d the cor- 
responding numerical results are shown in Fig. [2] 
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Figure 2. Relative finite volume corrections 
(Rm = — M v )/M v ) to the pion mass using 

the Luscher formula @ or the modified version 
which includes subleading exponentials I jlljl and 
the chiral expansion of the itit amplitude up to 
NNLO. 



A somewhat surprising feature of the numerical 
results is that the NLO corrections are quite sub- 
stantial and of the order of 50% of the LO ef- 
fect. The comparison of NLO and NNLO however 
shows that chiral expansion does have a good con- 
verging behaviour for all values of and L in 
the region of applicability of the formula. 

A comparison to the numerical values obtained 
with the full one-loop CHPT formula Eq. 
on the other hand shows that for not so large 
values of M n L exponentially subleading terms 
may be important. As was suggested in JTj 
the most reliable estimate of the finite volume 
effects for the pion mass is obtained by adding 
to the full one-loop CHPT result the NLO and 
NNLO corrections obtained through the Luscher 
formula. There is actually a simple extension of 
the Luscher formula that provides a convenient 
combination of both approaches. The modified 
formula is the following: 

fc=i v 

dyF{zy)e-^ M2 ^ L + 0( e - m ) 

J — oo 



and can be derived with the following reasoning. 
In Luscher's paper a substantial part of the 
proof of the formula is devoted to showing that 
the leading exponential corrections are obtained 
from loop graphs where only one pion propagator 
is taken in finite volume (i.e. is made periodic, 
as in Eq. ©) and all other propagators are the 
standard, infinite-volume ones. Since he concen- 
trates only on the leading exponential correction, 
he then immediately drops all terms with \n\ > 1 
in the infinite sum which makes the propagator 
periodic. But nothing forbids to keep all terms in 
the infinite sum Eq. JJjJ and to follow the rest of 
the derivation which leads to Eq. (|I1|) . 
A few comments are in order: 

1. although Eq. Ijlll) contains exponentially 
subleading terms, it is not algebraically bet- 
ter than Eq. © because both give only the 
leading term exactly - numerically, however, 
the effect of subleading terms is important 
for moderate values of M n L as shown in 

Fig.m 

2. an inclusion of exponentially subleading 
terms which is similar in spirit to the one 
proposed here has been already introduced 
by Luscher by extending the integration 
boundaries to infinity in Eq. J^J; 

3. at the two- loop level there will be two types 
of contributions which are not included in 
Eq. (|11() : contributions where two pion 
propagators are taken in finite volume, and 
contributions to the loop integration from 
other singularities other than the pion pole 
in the propagator (note that both contribu- 
tions are absent at the one-loop level: in- 
deed if one inserts the leading order ampli- 
tude HI in Eq. one recovers exactly 
the full one-loop result 0). Only a com- 
plete two-loop calculation will show how big 
these contributions are numerically |29|. 

3.3. Decay constants in finite volume 

The discussion in the previous subsection has 
shown the usefulness of the asymptotic formula 
of Luscher in evaluating finite- volume corrections 
for masses, and calls for extensions to other quan- 
tities. Recently, an extension to decay constants 
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Figure 3. Relative finite volume corrections to the 
pion decay constant (Rf„ = {Ftt,l — F n )/F n ). 



has been provided in 13 and reads 



/oo 
dye-V M '+v 2L N F (iy) + 0(e- ML ) , (12) 
-oo 

where Np(i>) is (as one intuitively expects) re- 
lated to the (0|v4 M |37r) amplitude. In this case 
there is a subtlety related to the fact that the 
latter amplitude has a pole due to the direct cou- 
pling of the axial current to a pion, which then 
rescatters into three pions. This pole appears ex- 
actly in the kinematical region where it is needed 
in Eq. I|12|) and must be subtracted. The pre- 
scription for the subtraction of this pole has been 
discussed in detail in J3j , together with the phys- 
ical interpretation of the subtraction. 

This asymptotic formula can now be used for a 
numerical analysis in the same way as the Luscher 
formula for the pion mass: in the present case the 
needed infinite-volume amplitude is the one rel- 
evant for the decay of the r into a neutrino and 
three pions and has been calculated to one loop in 
CHPT |30|. The numerical results, including the 
effects of subleading exponential terms, included 
according to the analogous of Eq. (fTT|> are shown 
in Fig. |21 Like Liischer's formula, this extension 
to decay constants can also be applied to heavier 
hadrons: what determines the finite- volume cor- 
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Figure 4. Relative finite volume corrections to the 
kaon decay constant (Rf k = (Fk,l — Fk)/Fk)- 



rection to the decay constant of the hadron H is 
the amplitude (Q\A^\2ttH) . For the case of the 
kaon, e.g., this is the amplitude relevant for K14 
decays, which is nowadays known to two loops 
|31|. A numerical analysis based on the one-loop 
evaluation of this amplitude gives the results 
shown in Fig. Finite volume corrections for 
F„ and Fk are also of phenomenological interest, 
because the lattice calculation of these two quan- 
tities can be used to determine V us , as suggested 
by Marciano [3333 1 an d performed by the MILC 
collaboration [31;. A complete numerical analy- 
sis of the masses and decay constants of all the 
octet of pseudoscalars is in progress |29| . 

4. Conclusions 

There has been a lot of activity recently in cal- 
culating finite volume effects analytically. I have 
briefly reviewed the theoretical tools needed for 
performing these calculations as well as some of 
this recent activity. The cost of a lattice calcula- 
tion grows very rapidly with the box size - the cal- 
culations discussed here show that in many cases 
it is unnecessary to make the infinite-volume ex- 
trapolation numerically, and one can correct for 
finite- volume effects analytically. 
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